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Abstract 

We consider the coupled Einstein-Dirac-Maxwell equations for a static, spherically 
symmetric system of two fermions in a singlet spinor state. Soliton-like solutions are 
constructed numerically. The stability and the properties of the ground state solutions 
are discussed for different values of the electromagnetic coupling constant. We find 
solutions even when the electromagnetic coupling is so strong that the total interaction 
is repulsive in the Newtonian limit. Our solutions are regular and well-behaved; this 
shows that the combined electromagnetic and gravitational self-interaction of the Dirac 
particles is finite. 

1 Introduction 

In the recent paper we studied the coupled Einstein-Dirac (ED) equations for a spheri- 
cally symmetric system of two fermions in a singlet spinor state. Using numerical methods, 
we found particle- like solutions and analyzed their properties. In this letter, we extend 
our results to the more physically relevant system where the Einstein-Dirac equations are 
coupled to an electromagnetic field. We find particle-like solutions which are linearly sta- 
ble (with respect to spherically symmetric, time-dependent perturbations), and have other 
interesting properties. 

The general Einstein-Dirac-Maxwell (EDM) equations for a system of n Dirac particles 

are 

i n 

R)--R8) = -8irT; , (G-m)* a = , V k F> k = 47re^^ a OT a , (1.1) 

1 o=l 

where Tj is the sum of the energy-momentum tensor of the Dirac particles and the 
Maxwell stress-energy tensor. Here the G 3 are the Dirac matrices, which are related 
to the Lorentzian metric via the anti-commutation relations 

gf( X ) = \ {&{X), G k (x)} , 

Fjk is the electromagnetic field tensor, and ^> a are the wave functions of fermions of mass m 
and charge e. The Dirac operator is denoted by G; it depends on both the gravitational and 
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the electromagnetic field (for details, see e.g. 0). In order to get a spherically symmetric 
system, we consider, as in two fermions with opposite spin, i.e. a singlet spinor state. 
Using the ansatz in Q, we reduce the Dirac 4-spinors to a 2-component real spinor system 
(a,f3). We show numerically that the EDM equations also have particle-like solutions, 
which are characterized by the "rotation number" n = 0,1,... of the vector (a, (5). In 
contrast to Q, we restrict our attention to studying the case n = 0, the ground state; this 
case illustrates quite nicely the physical effects due to the addition of the electromagnetic 
interaction. We anticipate that the situation for solutions with higher rotation number 
will be qualitatively similar. 

The relative coupling of the electromagnetic and the gravitational field is described 
by the parameter (e/m) 2 . If we consider the ground state solutions for fixed value of 
(e/m) 2 , we find that the mass-energy spectrum (i.e., the plot of the binding energy of the 
fermions vs. the rest mass) is a spiral curve which tends to a limiting configuration V. 
This implies the interesting result that for parameter values on T, there exist an infinite 
number of (n = 0)-solutions (the one of lowest energy is the ground state), while for 
parameter values near T, there are a large, but finite number of such solutions. For small 
coupling (i.e., for small m and (e/m) 2 < 1), our solutions are linearly stable with respect to 
spherically symmetric perturbations. If we compare the (ADM) mass p with the total rest 
mass 2m, we find that for small m, the total binding energy p — 2m is negative, implying 
that energy is gained in forming the singlet state. It is thus physically reasonable that 
such states should be stable. However, the stable solutions become unstable as the binding 
energy of the fermions increases; this is shown using Conley index methods together with 
bifurcation theory (see || Part IV]). 

In order to study the effect of the electromagnetic interaction in more detail, we look at 
the behavior of the solutions as the parameter (e/m) 2 is varied. For weak electromagnetic 
coupling (e/m) 2 <C 1, the solutions are well-behaved and look similar to the solutions of 
the ED equations For (e/m) 2 > 1, the form of the solutions changes drastically. In 
a simplified argument, this can already be understood from the nonrelativistic, classical 
limit of the EDM equations. Namely, according to Newton's and Coulomb's laws, the 
force between two charged, massive point particles has the well-known form 



m 2 e 2 . 
■— + - ( • ) 



(we work in standard units h = c = Q = 1). For (e/m) 2 < 1, the gravitational attraction 
dominates the electromagnetic repulsion, and it therefore seems reasonable that we get 
bound states. For (e/m) 2 > 1, however, the total force is repulsive, and clasically one 
can no longer expect bound states. The EDM equations, however, do have solutions even 
for (e/m) 2 > 1. This is a surprising effect which can again only be explained by the 
nonlinearity of Einstein's equations. For such solutions to exist, however, the rest mass m 
of the fermions must be sufficiently large. 

We remark that there are related works H p, M], where the authors obtain soliton 
solutions for the Dirac-Maxwell equations (in the absence of gravity), but these solutions, 
unfortunately, have the undesirable feature of having negative energy. 
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2 The Equations 

We choose polar coordinates (t, r, i?, <p) and write the metric in the form 

ds 2 = _L dt 2 _ _L dr 2 _ r 2 M 2 _ r 2 sin 2# dip 2 

with positive functions A(r), T(r). Using the ansatz from [jl|, Eqns. (3. 4), (3. 6)], we describe 
the Dirac spinors with two real functions a, (3. For the derivation of the corresponding 
EDM equations, we simply modify the ED equations jl], Eqns. (5.4)-(5.8)]. In the Dirac 
equation, the electromagnetic field is introduced by the minimal coupling procedure dj — > 
dj — ieAj, where e is the unit charge and A the electromagnetic potential (see In 
the static case, the fermions only generate an electric field; thus we can assume that the 
electromagnetic potential has the form A = (—</>, 0) with the Coulomb potential <fi. Since 
the time-dependence of the wave functions is a plane wave exp(— icut), minimal coupling 
reduces to the replacement u — > u — ecft. Thus the Dirac equations are 

VAa = - a - ((u - e<t>)T + m) p (2.1) 
r 

VAp' = ((w - e4>)T - m) a . (2.2) 



r 

In the Einstein equations, we must include the stress-energy tensor of the electric field Q; 
this gives 

r A' = 1 - A - 16vr(c; - e(f>)T 2 {a 2 + p 2 ) - r 2 AT 2 ((f>') 2 (2.3) 

T' 1 
2rA— = A - 1 - 16tt(lo - ec/))T 2 (a 2 + P 2 ) + 32tt - T aP + 16tt mT (a 2 - P 2 ) 
T r 

+r 2 AT 2 {<j)') 2 . (2.4) 

Using current conservation, Maxwell's equations VkF^ k = 4-7re J2a=l ^aG^^a reduce to 
the single second-order differential equation 



r 



r 
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r 2 A 4>" = -8vr e {a 2 + p 2 ) - \^2rA + r 2 A^+'— A' j $ . (2.5) 

The normalization condition for the wave functions Eqn. (5.8)] remains unchanged, 
namely 

oo r p 1 

{a 2 + p2 ) dr = (26) 

V A 47T 

We seek smooth solutions of these equations, which are asymptotically Minkowskian and 
have finite (ADM) mass p, 

lim T(r) = 1 (2.7) 

I — >oo 

p := lim - (1 - A(r)) < oo . (2.8) 

r — >oo 2 

Furthermore, we demand that the electromagnetic potential vanishes at infinity, 

lim Sir) = . (2.9) 

1 — *oo 

The equations ( |2.1[ )-( ^5| ) are invariant under the gauge transformations 

<f){r) -> 4>(r) + k , lo uj + en ,kEIR. (2.10) 
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3 Construction of the Solutions 



The construction of solutions is, (analogous to ]l|]), simplified by the following scaling 
argument. We weaken the conditions fl2.6|), fl2.7j), (P^S|) to 



^ / (a 2 +(3 2 ) —=dr < oo , ^ lim T(r) < oo , lim Mr) < oo , (3.1) 
Jo \/ A r ^°° r ^°° 

and set instead 

T(0) = 1 , 0(0) = , m = 1 . (3.2) 

This simplifies the discussion of the equations near r = 0; indeed, a Taylor expansion 
around the origin gives 

Q (r) = ai r + C(r 2 ) , /3(r) = C(r 2 ) 

A(r) = 1 + 0(r 2 ) , T(r) = 1 + 0(r 2 ) 

(j){r) = 0{r 2 ) . 

The solutions are now determined by only three real parameters e, u>, and ct\. For a given 
value of these parameters, we can construct initial data at r = and, using the standard 
Mathematica ODE solver, we shoot for numerical solutions of the modified system fl2.l|) - 
( |2.5D , (3.2). By varying ui (for fixed e and ai), we can arrange that the spinors (a, (3) 
tend to the origin for large r. As one sees numerically, the so-obtained solutions satisfy 
the conditions ( fO| ) and ( |37T]) . 

For a given solution (a, f3, A, T, </>) of this modified system, we introduce the scaled 
functions 

a(r) = ^a(Ar) , (3(r) = ^/3(\r) 

A(r) = A{Xr) , f(r) = r™ 1 T(Ar) 

0(r) = r (p(Xr) . 

As one verifies by direct computation, these functions satisfy the original equations ( |2.lD - 
([2.8D if the physical parameters are transformed according to 



rh = Am , Co = At lo , e = Ae , 

where the scaling factors A and r are given by 

i 

A = (471-/ (a 2 + /3 2 ) -= dr , r = lim T(r) 

\ Jo x/A J r ^°° 



Then the condition ( |2.9|) can be fulfilled by a suitable gauge transformation ( p. 10 ). No- 
tice that (e/m) 2 = e 2 is invariant under the scaling. Therefore it is convenient to take 
(e/rh) 2 (and not e 2 itself) as the parameter to describe the strength of the electromagnetic 
coupling. 

We point out that this scaling technique is merely used to simplify the numerics; for 
the physical interpretation, however, one must always work with the scaled tilde solutions. 
Since the transformation from the un-tilde to the tilde variables is one-to-one, our scaling 
method yields all the solutions of the original system. We will from now on consider only 
the scaled solutions; for simplicity in notation, the tilde will be omitted. 
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Figure 1: Solution of the EDM equations for parameter values (e/m) 2 
0.7639, u = 0.6807, p = 1.15416 (a'(0) = 0.05361). 



0.7162, m 



4 Properties of the Solutions 

The solutions we found have different rotation number n = 0, 1, . . . of the vector (a,/3). 
In the nonrelativistic limit, n coincides with the number of zeros of the corresponding 
Schrodinger wave functions, and thus n = corresponds to the ground state, n = 1 to 
the first excited state,. . . , etc. Because of the nonlinearity of our equations, n does not in 
general have this simple interpretation. In the following, we will restrict ourselves to the 
n = solutions. A plot of a typical solution is shown in Figure [l|. 

For all considered solutions, the spinors (a, (3) decay exponentially at infinity. This 
means physically that the fermions have a high probability of being confined to a neigh- 
borhood of the origin. Since the spinors decay so rapidly at infinity, our solutions asymp- 
totically go over into spherically symmetric solutions of the Einstein-Maxwell equations; 
i.e. the Reissner-Nordstrom solution B. More precisely, the behavior for large r is 



A(r) ps T(r 



-2 



1 



2p (2eY 

' 2 



2c 
r 



Thus asymptotically, our solution looks like the gravitational and electrostatic field gen- 
erated by a point particle at the origin having mass p and charge 2e . In contrast to the 
Reissner-Nordstrom solution, however, our solutions have no event horizons or singulari- 
ties. This can be understood from the fact that we consider quantum mechanical particles 
(instead of point particles) , which implies that the wave functions are delocalized accord- 
ing to the Heisenberg Uncertainty Principle. As a consequence, the distribution of matter 
and charge are also delocalized, and this prevents the metric from forming singularities. 

The situation near the origin r = 0, on the other hand, is parametrized by the rest 
mass m and the energy uj of the fermions. In Figure [| the binding energy m — u is plotted 
versus m for different values of the parameter (e/m) 2 . One sees that m — to is always 
positive, which means that the fermions are in a bound state. For weak electromagnetic 
coupling (see Figure |2|, plots A and B), the curve has the form of a spiral which starts at 
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Figure 2: Binding Energy m — uj of the Fermions for (e/m) 2 = (A), 0.7162 (B), 0.9748 
(C), 1 (D), and 1.0313 (E). 



the origin. The binding energy becomes smaller for fixed m and increasing (e/m) 2 ; this 
is because the electromagnetic repulsion weakens the binding. The mass-energy spectrum 
when (e/m) 2 <C 1 has a similar shape as in the case without the electromagnetic interaction 
0. The stability techniques and results of [ffl] can be generalized directly: For small m, 
one can use linear perturbation theory to show numerically that the solutions are stable 
(with respect to spherically symmetric perturbations). The stability for larger values of 
m can be analyzed with Conley index theory (see ||), where we take m as the bifurcation 
parameter. The Conley index theory yields that the stability /instability of a solution 
remains unchanged if the parameter m is continuously varied and no bifurcations take 
place. Moreover, at the bifurcation points, the Conley index allows us to analyze the 
change of stability with powerful topological methods. We find that all the solutions on 
the "lower branch" of the spiral (i.e., on the curve from the origin up to the maximal value 
of m) are stable, whereas all the solutions on the "upper branch" are unstable. 

The form of the energy spectrum changes when (e/m) 2 ~ 1. This is the regime where 
the electrostatic and gravitational forces balance each other in the classical limit (|1.2j ). 
Since the assumption of classical point particles does not seem to be appropriate for our 
system, we replace ( |1.2|) by taking the nonrelativistic limit of the EDM equations more 
carefully. For this, we assume that m and e are small (for fixed (e/m) 2 ). The coupling 
becomes weak in this limit; thus A, T ~ 1, (f> 0. The Dirac equations give uj ~ m and 
a 3> (3. Therefore the EDM equations go over into the Schrodinger equation with the 
Newtonian and Coulomb potentials, 

f — — A + ecb + mV) ^ = E * 
V 2m / 

-AV = -8vrm|*| 2 , -A0 = 8vre|^| 2 

(with E = uj — m, 'l'(r) = a(r)/r, V = 1 — T; A equals the radial Laplacian in IR 3 ). One 
sees from these equations that, as in ( |1.2[ ), the Newtonian and Coulomb potentials are 
just multiples of each other; namely V = —m(/)/e. For (e/m) 2 > 1, the total interaction 
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Figure 3: Total Binding Energy p- 2m for (e/m) 2 = (A), 0.7162 (B), 0.9748 (C), 1 (D), 
and 1.0313 (E). 

is repulsive, and the Schrodinger equation has no bound states. We conclude that, in the 
limit of small to, the EDM equations cannot have particle-like solutions if (e/m) 2 > 1. In 
other words, the mass-energy plots of Figure ^ can only start at to = if (e/m) 2 < 1. 
This is confirmed by the numerics (see Figure plots C, D, and E). For (e/m) 2 = 1, the 
plot tends asymptotically to the point (to = do, to — uj = 0). It is surprising that we still 
have bound states for (e/m) 2 > 1. In this case however, there are only solutions if to is 
sufficiently large and smaller than some threshold value where the binding energy of the 
fermions becomes zero. 

In Figure ||, the total binding energy p — 2m is plotted for different values of (e/m) 2 . 
For (e/m) 2 < 1, p — 2m is negative for the stable solutions. For (e/m) 2 > 1, however, 
p — 2m is always positive. This indicates that these solutions should be unstable, since 
one gains energy by breaking up the binding. 
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